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Abstract  Coarse-grained  molecular  statics/dynamics 
methods  are  first  used  to  investigate  degradation  in  the 
PPTA  fiber/yam  tensile  strength,  as  a  result  of  the  prior 
compressive  or  tensile  loading.  PPTA  fibers/yarns  experi¬ 
ence  this  type  of  loading  in  the  course  of  a  plain- weaving 
process,  the  process  which  is  used  in  the  fabrication  of 
ballistic  fabric  and  flexible  armor.  The  more  common  all¬ 
atom  molecular  simulations  were  not  used  to  assess 
strength  degradation  for  two  reasons:  (a)  the  size  of  the 
associated  computational  domain  rendering  reasonable 
run-times  would  be  too  small  and  (b)  modeling  of  the 
mechanical  response  of  multi-fibril  PPTA  fibers  could  not 
be  carried  out  (again  due  to  the  limited  size  of  the  com¬ 
putational  domain).  However,  all-atom  simulations  were 
used  to  (a)  define  the  coarse-grained  particles  (referred  to 
as  “beads”)  and  (b)  parameterize  various  components  of 
the  bead/bead  force-field  functions.  In  the  second  portion 
of  the  work,  a  simplified  finite-element  analysis  of  the 
plain-weaving  process  is  carried  out  in  order  to  assess 
the  extent  of  tensile- strength  degradation  in  warp  and 
weft  yarns  during  the  weaving  process.  In  this  analysis,  a 
new  material  model  is  used  for  the  PPTA  fibers/yarns. 
Specifically,  PPTA  is  considered  to  be  a  linearly  elastic, 
transversely  isotropic  material  with  degradable  longitudi¬ 
nal-tensile  strength  and  the  longitudinal  Young’s  modulus. 
Equations  governing  damage  and  strength/stiffness 
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degradation  in  this  material  model  are  derived  and 
parameterized  using  the  coarse-grained  simulation  results. 
Lastly,  the  finite-element  results  are  compared  with  their 
experimental  counterparts,  yielding  a  decent  agreement. 

Introduction 

The  present  work  deals  with  the  damage  induced  into 
/^-phenylene  terephthalamide  (PPTA)  fibers/yarns  during 
the  plain- weaving  process,  the  process  which  is  used  in  the 
fabrication  of  flexible  armor.  Since  the  plain-weaving 
process  is  usually  based  on  the  use  of  two  orthogonal 
families  of  yarns  (generally  referred  to  as  warp,  or  end,  and 
weft,  fill  or  pick  yams),  the  work  addresses  differences  in 
the  mechanical  history  of  the  two  yarn  families  and  the  role 
of  these  differences  in  the  extent  of  damage  sustained  by 
the  two  yarn  families.  Consequently,  the  main  topics 
overviewed  in  the  remainder  of  this  section  include  (a)  the 
basics  of  plain- weaving;  (b)  degradation  of  weft  and  warp 
yarns  during  weaving;  (c)  PPTA  high-performance  fibers/ 
yarns;  and  (d)  prior  experimental  investigations  of  weav¬ 
ing-induced  yarn  damage. 

The  basics  of  plain-weaving 

Figures  la-d  show  the  most  common  (plain,  basket,  satin, 
and  triaxial)  weaving  patterns  used  in  flexible-armor 
applications.  Plain- weaving.  Fig.  la,  is  the  subject  of  the 
present  work.  Plain-weaving  is  one  of  the  key  fabric-pro¬ 
duction  technologies.  It  involves  interlacing  two  families  of 
yarns/fibers,  in  a  predefined  order,  so  as  to  form  a  fabric  of 
the  desired  architecture  and  properties.  Weaving  is  carried 
out  using  a  machine  called  a  loom.  The  warp  yarns,  which 
form  the  basic  structure  of  the  fabric  and  are  aligned  in  the 
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Fig.  1  Most  common  weaving  patterns  used  in  flexible-armor  applications:  a  plain;  b  basket;  c  satin;  and  d  triaxial 


lengthwise  direction  of  the  fabric,  run  from  the  front  to  the 
back  of  the  loom.  On  the  other  hand,  the  weft  yarns  run 
orthogonally  to  the  warp  yarns,  and  the  two  sets  of  yams 
are  interlaced. 

Typically,  the  weaving  process  involves  repeated 
application  of  the  following  four  basic  steps: 

(a)  shedding — within  this  step,  warp  yarns  are  separated 
into  two  groups/sheets,  with  the  pattern  of  separation 
depending  on  the  weave  architecture,  e.g.,  for  a 
plain-weave,  alternating  warp  yarns  are  placed  into 
the  same  sheet.  Separation  between  the  two  sheets  is 
referred  to  as  shed,  and  its  size  has  to  be  sufficiently 
large  to  enable  the  passage  of  the  weft-yarn  carrier/ 
shuttle; 

(b)  picking — within  this  step,  weft/pick  yams  are  laid 
between  the  two  warp  sheets,  i.e.,  into  the  shed; 

(c)  heating  in — within  this  step,  a  comb-like  tool  called 
a  reed  is  used  to  push  weft  yarns  against  the  fabric 
already  formed;  and 

(d)  taking  up  and  letting  off — within  this  step,  the  fabric 
is  being  wound  onto  a  cloth-/take-up-roller,  while 
warp  yarns  are  being  'det  in''  from  the  weaver’s 


beam.  Strictly  speaking,  this  is  not  a  separate  step, 
but  a  step  that  is  taking  place  at  a  constant  rate. 

Once  one  four-step  cycle  is  completed,  the  action  is 
returned  to  the  shedding  step,  within  which  a  new  shed  is 
formed  in  accordance  with  the  desired  weave  architecture.  In 
other  words,  some  or  all  of  the  warp  yams  previously  residing 
in  one  sheet  are  moved  to  the  previous  position  of  the  other 
sheet.  In  the  case  of  plain-weave,  the  entire  two  warp-yam 
sheets  exchange  their  places.  In  this  process,  the  weft-yarn  is 
sandwiched  between  two  sets  of  crimped/undulated  warp 
yams,  and  can  experience  minor  crimping/waviness.  This 
description  clearly  indicates  that  warp-yams  experience,  in 
general,  more  deformation/crimping  (and,  hence,  potentially 
more  damage)  during  the  weaving  process,  than  their  weft 
counterparts.  Consequently,  when  the  fabrics  are  made  of 
ballistic-class  yams,  the  warp  yams,  in  the  woven  fabric, 
generally  possess  inferior  mechanical  properties. 

Degradation  of  weft  and  warp  yarns  during  weaving 

The  analysis  presented  in  the  previous  section  clearly 
revealed  that,  during  the  weaving  process,  yarns  (in 
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particular,  the  warp  yarns)  can  be  subjected  to  a  consid¬ 
erable  amount  of  crimping.  Since  crimping  is  a  bending 
process,  the  largest  strains/stresses  associated  with  it  reside 
on  the  surface  of  the  yarns.  When  these  strains/stresses  are 
sufficiently  high,  surface  damage  can  be  introduced  into  the 
yarns,  causing  degradation  in  their  mechanical  properties. 

An  additional  mode  of  surface  damage  arises  from  the 
fact  that  during  the  weaving  process,  crossing  yarns  in 
contact  undergo  relative  sliding.  This  process  can  cause 
surface  abrasion  in  the  contacting  yarns,  another  mode  of 
yarn- surface  damage.  This  damage  mode  can  further 
compromise  the  structural  integrity  and  mechanical  prop¬ 
erties  of  the  yarns  residing  in  the  woven  fabric. 

Also,  before  being  used  in  the  weaving  process,  yams 
are  often  coated  for  surface  protection  and,  hence,  the 
removal  of  the  coating  from  the  woven  fabric  requires 
fabric  washing  and  drying.  These  processes  can  also  impart 
some  damage  (not  only  of  a  surface-  but  also  of  a  bulk- 
type)  within  the  yarns  residing  in  the  woven  fabric. 

Based  on  the  analysis  presented  above,  the  following 
main  observations  can  be  made:  (a)  yams  residing  in  the 
woven  fabric  possess  mechanical  properties  which  are 
inferior  relative  to  their  counterparts  in  the  virgin  yarn  and 
(b)  the  degree  of  this  inferiority  is  expected  to  be  larger  in 
the  warp  yarns. 

PPTA  high-performance  fibers/yams 

Since  the  formulation  and  construction  of  the  molecular- 
level  PPTA  computational  models  used  in  the  present  work 
entailed  detailed  knowledge  of  the  basic  PPTA  chemistry, 
bonding,  crystal,  and  hierarchical  structures  and  fabrica¬ 
tion-induced  defects,  a  brief  overview  of  these  material 
aspects  is  presented  in  the  remainder  of  this  section. 

Class  identification 

PPTA  materials  (e.g.,  Kevlar®,  Twaron®’  Technora®  etc.) 
fall  into  the  category  of  aromatic  polyamide  or  aramid 
materials. 

PPTA  molecule s/chains 

The  PPTA  repeat  units/mers  appears  in  two  distinct: 
(a)  trans  {''on  the  opposite  side'')  and  (b)  cis  {"on  the  same 
side")  molecular-level  conformations.  These  conforma¬ 
tions  are  shown  using  a  ball-and-stick  representation  in 
Fig.  2a-b,  respectively.  To  help  with  interpretation  of  the 
molecular  stmcture,  the  chemical  elements  are  labeled  in 
Fig.  2a.  Examination  of  Fig.  2a-b  reveals  that  the  repeat 
unit  consists  of  two  phenylene  rings/moieties  (C6H4RR’) 
connected  by  an  amide  linkage  (-CO-NH-).  Although 
the  PPTA  molecule  can  assume  either  the  trans  or  cis 


stereo-isomeric  conformations,  the  former  conformation  is 
almost  exclusively  observed.  The  reason  for  this  predom¬ 
inance  is  that  the  trans  conformation  promotes  the 
stretching/extension  of  PPTA  molecules,  lowering  the 
system  energy.  In  addition,  extended  molecules  can  readily 
align,  forming  sheets  consisting  of  nearly  parallel,  regu¬ 
larly  spaced  PPTA  chains.  Formation  of  these  sheets  is 
associated  with  further  reductions  in  the  system  energy. 
Regular  stacking  of  these  sheets  leads  to  the  formation  of 
low-energy  crystalline  PPTA  fibrils  (subfiber  units,  elon¬ 
gated,  and  coaxial  with  the  fibers).  In  sharp  contrast,  the  cis 
conformation  hampers  the  formation  of  these  low-energy 
structures,  which  explains  its  rare  occurrence. 

PPTA  sheets 

Due  to  a  large  difference  in  electronegativity  between 
oxygen  and  hydrogen,  the  electron  charge  density  within 
the  PPTA  unit  cell  is  increased  around  the  oxygen  atoms 
and  lowered  around  the  hydrogen  atoms,  giving  rise  to  the 
formation  of  (relatively  strong)  dipoles.  The  presence  of 
dipoles  along  with  (i)  the  lone  electron  pairs  on  the  oxygen 
atoms;  (ii)  close  proximity  of  the  oxygen  and  hydrogen 
atoms  of  neighboring  parallel  PPTA  molecules  within  the 
same  sheet;  and  (iii)  the  operation  of  London-dispersion 
forces  typically  results  in  the  formation  of  (strong) 
hydrogen  bonds.  London-dispersion  forces  are  temporary 
intermolecular  forces  which  are  caused  by  instantaneous 
polarization  of  the  interacting  molecules  due  to  repulsion 
of  their  adjacent  electrons  [1].  Within  the  hydrogen  bonds, 
the  N-H  species  acts  as  a  bond  donor,  while  O  acts  as  a 
bond  acceptor.  The  high  strength  of  hydrogen  bonds  stems 
from  the  covalent-like  bonding  between  the  donor  and 
acceptor  atoms  (accompanying  the  redistribution  of  the 
lone  electron-pair  charge  density),  electrostatic  attractions 
between  the  opposite  poles  of  the  neighboring  chains,  and 
the  development  of  London-dispersion  forces.  The  forma¬ 
tion  of  the  aforementioned  sheet-like  structures  is  believed 
to  be  the  result  of  strong  lateral  hydrogen  bonding  between 
adjacent  parallel  PPTA  molecules. 

PPTA  fibrils 

As  mentioned  above,  stacking  of  sheets  results  in  the  for¬ 
mation  of  crystalline  fibrils,  the  axes  of  which  are  aligned 
with  the  fiber  axis.  However,  orientation  of  the  fibril  crystal 
structure  within  a  plane  normal  to  the  fiber  axis  varies 
among  the  fibrils.  In  other  words,  PPTA  fibrils  within  a 
single  fiber  resemble  columnar  grains  observed  in  direc¬ 
tionally  solidified  as-cast  metallic  micro  structures.  How¬ 
ever,  in  contrast  to  the  columnar  metallic  grains,  the  PPTA 
fibrils  are  generally  found  to  possess  a  superstructure 
consisting  of  small  amplitude  (ca.  20  nm),  high  wavelength 
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Fig.  2  a  Trans--,  and  b  cis- 
molecular  conformations  in 
typical  PPTA-based  polymeric- 
material  molecules/chains; 
c  NH2-type  chain-end;  and 
d  S03H-type  side-group  defects 
present  in  PPTA  fibrils/fibers 


chain 

defect 


SO3H  side-group 
defect 


(ca.  300-500  nm)  accordion-style  ‘‘pleats.  ”  The  presence 
of  these  pleats  is  often  cited  as  one  of  the  causes  for  the 
superior  tensile  ductility  of  PPTA  fibrils/fibers  [e.g.,  2]. 
While  intra- sheet  bonding  is  of  a  strong  hydrogen-bonding 
type,  inter- sheet  bonding  involves  weak  van  der  Waals 
interactions  and  the  formation  of  Ti-type  molecular  orbitals. 
Thus,  PPTA  fibrils  possess  a  substantial  level  of  transverse 
anisotropy. 

In  addition  to  the  stimulating  effect  of  strong  hydrogen 
bonding,  PPTA  sheet  and  fibril  formation  are  also  believed  to 
be  influenced  by  the  unique  molecular  (or,  more  specifically, 
back  bone)  structure  of  PPTA  chains,  the  structure  which  is 
associated  with  a  high  value  of  the  bending  rigidity/stiffness. 
Specifically,  due  to  limited  molecular  flexibility  and  the 


absence  of  molecular  folding  of  PPTA  chains,  the  PPTA 
fibril/fiber  micro  structure  is  not  the  commonly  observed 
mixture  of  crystalline  and  amorphous  phases.  Rather,  the 
microstructure  of  PPTA  fibrils/fibers  is  characterized  as 
being  either  paracrystalline  or  fully  crystalline.  While  in  both 
structures,  molecules  are  aligned  with  the  fiber  axis,  the 
sheets  are  stacked  with  a  significantly  higher  degree  of  reg¬ 
ularity  in  the  fully  crystalline  than  in  the  paracrystalline 
structure.  Another  factor  often  cited  as  the  promoter  of  the 
fully  crystalline/paracrystalline  structure  within  the  PPTA 
fibrils/fibers  is  the  presence  of  the  planar  phenylene  and 
amide  groups  within  the  PPTA  chains. 

As  mentioned  above,  inter-sheet  bonding  within  PPTA 
fibrils  is  relatively  weak  and,  consequently,  fibrils  tend  to 
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Fig.  3  Topological  relationships  between:  a  fibrils  within  a  single 
fiber;  and  b  fibers  within  a  lightly  twisted  yam 


possess  various  crystallographic  and  topological  inter-sheet 
defects  (e.g.,  local  delaminations  13]).  This  is  the  reason 
that  PPTA  fibrils  tend  to  buckle  during  axial  and  lateral 
(sheet  edge-on)  compressive  loading,  and  to  possess  rela¬ 
tively  inferior  axial/transverse  compressive  properties. 
Weak  inter-sheet  bonding  can  also  result  in  increased  water 
absorption,  further  compromising  structural  integrity  and 
durability  of  PPTA  fibrils  12]. 

PPTA  fibers 

As  mentioned  above,  PPTA  fibers  are  essentially  assem¬ 
blies  of  parallel  fibrils  (nearly  coaxial  with  the  fiber  itself). 
While  it  is  generally  believed  that  fibril  bundles  within  the 
fiber  are  held  together  by  weak  van  der  Waals  or  Coulomb 
forces,  it  is  possible  that  chain  cross-linking  and  the  pre¬ 
sence  of  tie-fibrils  can  play  a  role  in  the  formation  of  fibril 
bundles.  While  all  fibrils  contain  the  same  crystalline  or 
nearly  crystalline  structure,  they  possess  different  relative 
crystallographic  orientations.  A  schematic  of  the  basic 
architecture  of  the  single  PPTA  fiber  and  its  relationship  to 
the  constituent  fibrils  is  depicted  in  Fig.  3a.  Examination  of 
this  figure  reveals  that  the  fibril  bundles  form  the  fiber  core 
and  that  the  fibers  also  contain  a  relatively  thin  skin  region. 
While  this  region  is  comprised  of  highly  oriented  PPTA 
chains,  it  is  usually  of  a  non-crystalline  character.  It  should 
be  noted  that,  while  the  schematic  displayed  in  Fig.  3a 
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shows  only  six  fibrils,  actual  PPTA  fibers  typically  reveal 
the  presence  of  hundreds  of  such  fibrils  within  their  core  (at 
a  given  longitudinal  location). 

PPTA  yarns 

PPTA  yarns  are  considered  as  assemblies  of  nearly  parallel 
PPTA  fibers  which  are  mechanically  engaged  by  either  the 
application  of  a  light  twisting  to,  or  by  wrapping  a  fiber 
around,  the  fiber  assembly.  In  general,  neither  of  these  two 
ways  of  yam  fabrication  gives  rise  to  a  significant  damage 
to  the  fiber  assembly.  A  schematic  of  the  basic  architecture 
of  the  single  “twisted”  PPTA  yarn  and  its  relationship  to 
the  constituent  fibers  is  depicted  in  Fig.  3b. 

Role  of  defects/fiaws 

Since  PPTA  fibers  consist  of  bundles  of  nearly  parallel  and 
coaxial  (fully  crystalline  or  paracrystalline)  fibrils,  their 
properties  (as  is  the  case  for  most  crystalline  materials)  are 
greatly  affected  by  the  presence,  the  type,  and  the  con¬ 
centration  of  various  crystallographic/topological  defects. 
For  example,  within  the  crystalline  PPTA  fibrils,  sheet 
stacking  is  typically  of  an  ABABAB  type.  This  stacking  is 
disturbed  in  the  case  of  paracrystalline  fibrils,  leading  to 
the  formation  of  higher  energy  stacking  faults.  Crystal¬ 
lographic/topological  defects  do  not  only  affect  fiber 
properties  but  also  the  properties  and  the  performance  of 
the  coarser-scale  stmctures  (e.g.,  yarns,  fabrics,  plies, 
laminae,  and  laminates).  The  type,  size/potency,  and  the 
concentration  of  these  flaws  are  a  sensitive  function  of  the 
PPTA  synthesis  and  fiber  fabrication  processes.  A  sum¬ 
mary  of  the  processes  of  PPTA  synthesis  and  fiber  fabri¬ 
cation,  as  reported  in  the  open  literature  14-6],  was 
provided  in  our  recent  work  17-12].  The  summary  clearly 
established  the  effect  of  the  synthesis/processing  parame¬ 
ters  on  the  type,  size,  and  the  concentration  of  the  crys¬ 
talline  and  topological  defects  within  the  PPTA  fibers.  In 
the  present  work,  the  effect  of  two  types  of  PPTA  defects, 
the  so-called  “chain-end”  and  the  so-called  “side-group,” 
is  investigated.  Molecular-level  structures  of  two  small 
regions  within  the  PPTA  fibril,  each  containing  one  of 
these  defects  are  depicted  in  Fig.  2c-d.  In  the  case  of 
Fig.  2c,  the  chain-end  is  terminated  with  a  NH2  group, 
while  in  the  case  of  Fig.  2d,  the  inter- sheet  side-group 
consists  of  a  SO3H  molecule. 

Prior  experimental  investigations  of  weaving-induced 
yarn  damage 

Lim  et  al.  113]  investigated  the  damage  induced  into  the 
weft  and  warp  PPTA  yarns  during  weaving  and  the  bal¬ 
listic-penetration  resistance  of  two  PPTA  yarn-based 
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fabrics.  Their  results  clearly  revealed  that  (a)  the  warp 
yarns  extracted  from  the  fabric  possess  inferior  strength 
properties  to  their  weft  counterparts;  (b)  the  extent  of  warp- 
yarn  strength  degradation  increases  with  the  linear  density 
of  the  weft  yams.  This  observation  was  linked  to  the  fact 
that,  as  the  linear  density  of  the  weft  yarns  increases,  the 
extent  of  warp-yarn  crimping/waviness  (and,  hence,  dam¬ 
age)  also  increases;  and  (c)  the  ballistic  performance  of  the 
fabric  initially  increases,  and  peaks  then  begins  to  degrade 
as  the  weft-yarn  linear  density  increases.  These  findings 
have  been  rationalized  by  the  fact  that  initially  the  effect  of 
added  weft  yarns  has  a  dominant  effect,  while  at  higher 
values  of  the  weft-yarn  linear  density,  the  effect  of  warp- 
yarn  degradation  begins  to  prevail.  In  other  words,  there  is 
an  optimum  value  of  the  weft-yam  linear  density  associ¬ 
ated  with  the  maximum  ballistic-penetration  resistance  of 
the  fabric. 

Sanborn  et  al.  [14]  compared  the  mechanical  behavior  of 
warp  and  weft  yarns  extracted  from  plain-woven  Kevlar 
KM2  (a  commercial  grade  of  PPTA)  fabric  with  the  cor¬ 
responding  reference/non- woven  yarns.  Quasi-static  and 
intermediate  rate  experiments  were  conducted  using  a  Bose 
ElectroForce  Test  System,  while  high-rate  performance 
was  measured  using  a  miniature  tension  Kolsky  bar.  The 
results  obtained  revealed  that  (a)  warp  yarns  possess  infe¬ 
rior  strength  properties  with  respect  to  the  reference  yams 
and  (b)  the  weft  yarns,  in  a  statistical  sense,  possess  the 
same  strength  properties  as  the  reference  yams. 

Rocher  et  al.  [15]  investigated  the  extent  of  weaving- 
induced  strength  degradation  in  fabrics  made  of  hybrid 
polypropylene/glass  yarns,  at  different  testing  strain  rates. 
The  results  obtained  relatively  little  damage  (and,  hence, 
strength  degradation)  in  weft  yarns.  On  the  other  hand,  the 
fracture  strength  and  Young’s  modulus  of  the  warp  yams 
were  reduced,  on  average,  by  28  and  17  %,  respectively, 
relative  to  their  reference/non- woven  counterparts.  The 
work  also  showed  that  as  the  linear  density  of  the  weft 
yarns  increases  and,  hence,  the  extent  of  warp-yarn 
crimping  increases,  the  strength  and  stiffness  degradations 
in  the  warp  yarn  also  increase. 

Hunt  et  al.  [16,  17]  carried  out  an  extensive  investiga¬ 
tion  of  the  degradation  of  warp  and  weft  yams  during  the 
fabric-weaving  process,  as  well  as  the  degradation  of 
through-the-thickness  yarns  used  for  stitching  multi-ply 
fabric  structures.  A  variety  of  yams  were  examined, 
including  Dyneema  [highly  oriented  ultra-high  molecular 
weight  polyethylene  (UHMWPE)  filaments],  Vectran 
(polyester),  S-2  glass,  Kevlar  KM2,  etc.  The  results 
obtained  clearly  established  that  (a)  the  dominant  strength¬ 
degrading  mechanism  may  differ  from  one  yam  type  to  the 
other.  For  example,  in  Kevlar  KM2  and  Vectran,  the 
crimping-induced  damage,  in  the  form  of  kink  bands 
(primarily  in  the  warp  yams),  has  the  prevailing  effect. 


On  the  other  hand,  in  Dyneema,  surface  abrasion  appears  to 
be  the  dominant  fabric -fabrication  process-induced  yarn 
strength-degradation  mechanism;  (b)  in  multi-ply  fabric 
structures,  through-the-thickness  stitching  yarns  may 
experience  even  larger  strength  degradation  than  warp 
yarns  during  fabric- weaving;  and  (c)  among  the  yams 
investigated,  Dyneema  appears  to  have  the  highest  strength 
retention  when  used  as  cross-stitching  yarn. 

Main  objective 

The  two  main  objectives  of  the  present  work  are  (a)  to 
quantify  computationally  the  effect  of  prior  compressive 
and  tensile  (bending)  stresses  on  the  longitudinal-tensile 
strength  of  PPTA  fibers  using  coarse-grained  molecular 
methods  and  tools  and  (b)  to  assess  the  extent  of  strength 
degradation  experienced  by  the  warp  and  weft  yarns  during 
the  plain-weaving  process,  by  carrying  out  a  simplified 
finite-element  analysis  of  this  process. 

Paper  organization 

Details  regarding  the  coarse-grained  techniques  used 
including  the  construction  of  the  computational  model, 
formulation  and  parameterization  of  the  force-fields, 
computational  algorithms,  problem  formulation,  post-pro- 
cessing  data-reduction  analyses,  and  the  key  results 
obtained  are  all  presented  in  “Coarse-grained  fiber-scale 
computational  analysis”  section.  The  finite-element  ana¬ 
lysis  (including  the  results)  of  the  plain-weaving  process 
and  the  accompanying  degradation  in  the  warp-/weft-yarn 
strengths  is  described  in  “Finite-element  analysis  (FEA)  of 
plain- weaving  process”  section.  The  key  findings  resulting 
from  the  present  work  are  summarized  in  “Conclusions” 
section. 


Coarse-grained  fiber-scale  computational  analysis 

As  mentioned  earlier,  coarse-grained  fiber- scale  computa¬ 
tional  methods  have  been  employed  in  the  present  work  in 
order  to  investigate  degradation  in  the  PPTA  fiber/yarn 
tensile  strength,  as  a  result  of  the  prior  compressive  or 
tensile  loading.  PPTA  fibers/yams  experience  this  type  of 
loading  in  the  course  of  the  plain-weaving  fabric-manu¬ 
facturing  process. 

Within  the  coarse-grained  computational  framework, 
atoms,  ions,  and  inter-atomic  bonds  are  not  modeled 
explicitly  but  rather,  individual  atoms/ions  are  combined 
into  larger  single  particles,  referred  to  as  ''beads,''  where 
each  bead  represents  the  collective  degrees  of  freedom  of 
the  constituent  atoms/ions.  Neighboring  beads  residing 
within  the  same  chain  are  then  connected  using  bead/bead 
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bonds  (also  referred  to  as  connectors),  while  non-bonded 
beads  are  allowed  to  interact  via  van  der  Waals-like  and 
Coulomb-like  force-fields  (the  term  ''force-field''  will  be 
defined  later).  Since  the  coarse-grained  computational 
framework  employs  larger  particles  (relative  to  the  particle 
sizes  encountered  in  the  so-called  all-atom  computational 
analyses),  larger  computational  domains  can  be  analyzed 
(simulating  the  maximum  number  of  particles  permitted  by 
the  available  computational  and  storage  resources).  In 
addition,  by  integrating  out  the  high-frequency  vibrations 
of  the  constituent  atoms/ions  and  allowing  the  system 
dynamics  to  be  controlled  by  the  compliant/soft,  bead/bead 
interactions  (rather  than  by  stiff/hard,  atom/atom  interac¬ 
tions),  a  significantly  longer  time  scale  is  attained  (relative 
to  that  which  could  be  afforded  by  the  all-atom  computa¬ 
tional  method). 

It  should  be  noted  that,  while  only  coarse-grained 
molecular  methods  and  tools  are  used  in  the  present  work, 
the  results  of  the  all-atom  molecular-dynamics  simulations 
carried  out  in  our  recent  work  17-11]  were  utilized.  Spe¬ 
cifically,  all-atom  molecular-dynamics  calculations 
described  in  Refs.  17-11]  were  used  (a)  to  identify  the 
nature  and  the  number  of  types  of  coarse-grained  beads; 
(b)  for  each  bead  type,  to  quantify  the  bead  mass  and  bead 
size;  and  (c)  to  help  parameterize  intra-  and  inter-chain 
coarse-grained  force-fields  (a  set  of  mathematical  expres¬ 
sions  describing  various  contributions  to  the  potential 
energy  function  for  the  system  of  interacting  particles,  i.e., 
beads,  in  the  present  case).  A  brief  description  of  the  all¬ 
atom  computational  procedures  employed  in  our  recent 
work  17-11],  and  of  the  associated  results  used  in  the 
present  work,  is  provided  below. 

PPTA  coarse-graining  and  force-field  parameterization 

Within  the  all-atom  computational  methods,  molecular- 
level  models  of  individual  PPTA  chains  based  on  the  trans 
molecular  conformation  of  the  repeat  unit.  Fig.  2a  are 
constructed  using  the  Visualizer  118]  program  from  Ac- 
celrys.  Typically,  PPTA  chains  consisting  of  20  repeat 
units  and  having  the  molecular  weight  of  ca.  4760  Da  are 
used.  Atom/atom  bonding  and  non-bonding  interactions  are 
quantified  using  the  so-called  “COMPASS”  (Condensed- 
phase  Optimized  Molecular  Potentials  for  Atomistic  Sim¬ 
ulation  Studies)  force-field  119,  20].  This  highly  accurate 
force-field  is  of  an  ab  initio  type  since  most  of  its  param¬ 
eters  were  determined  by  matching  the  predictions  made  by 
the  ab  initio  quantum  mechanics  calculations  to  the  con¬ 
densed-matter  experimental  data.  The  force-field  defines 
mathematical  expressions  for  the  energy  contribution  of, 
and  the  force  acting  upon,  a  bead  as  a  function  of  its  local 
environment  (i.e.,  the  nature,  the  position,  and  the  bond 
status  of  the  neighboring  beads). 

^  Springer 


PPTA  coarse-graining 

Since  PPTA  can  be  considered  as  a  homopolymeric 
material,  and  in  order  to  be  able  to  handle  computational 
domains  of  a  larger  size,  only  one  type  of  bead  is  intro¬ 
duced  during  the  coarse-graining  procedure  of  the  PPTA 
chains.  Consequently,  the  bead  mass  was  set  equal  to  that 
of  the  PPTA  repeat  unit  while  the  volume  allotted  to  each 
bead  was  readily  derived  from  the  PPTA-fiber  mass 
density. 


Coarse-grained  force-field  functions 


The  coarse-grained  force-field  used  in  the  present  work 
includes  four  (functionally  distinct)  contributions  to  the 
system  potential  energy:  (a)  intra-chain  bond-stretch  term; 
(b)  inter-chain  two-bond-angle  term;  (c)  intra-sheet 
hydrogen-bond  term;  and  (d)  inter- sheet  van  der  Waals 
non-bond  interaction  term.  The  functional  forms  for  the 
four  potential-energy  terms  are  given  below: 

Intra-chain  bond-stretch  term  Covalent  bonding  within 
the  bead  type  PPTA  chains  is  represented  using  the  fol¬ 
lowing  form  of  a  harmonic- type  energy  function: 

1  ^ 

^stretch  _  _j^stretch^^  _  ^stretch^^ 


where  denotes  the  bond-stiffness  parameter,  r  the 

connector  (i.e.,  bead  bond)  length,  and  the  subscript  o  is 
used  to  denote  the  equilibrium  value  of  the  associated 
quantity.  Due  to  the  high  stretching  stiffness  of  the  PPTA 
molecules,  is  assigned  a  relatively  large  value,  as 

will  be  shown  later.  It  should  be  noted  that  the  bond- stretch 
term  describes  interactions  between  two  covalently  bonded 
beads  using  a  linear  axial  spring  model. 

Intra-chain  two -bond- angle  term  An  analogous  har¬ 
monic  function  is  used  to  represent  the  potential  energy 
contribution  associated  with  the  intra-chain  two-bond- 
angle  term  as  follows: 


jangle  _  _  jangle 


(2) 


where  denotes  the  bond-stiffness  parameter  and  9  the 
intra-chain  inter-bond  angle.  Due  to  the  high  bending 
rigidity  of  the  PPTA  molecules,  is  set  to  a  value  of  n, 
while  as  will  be  shown  later,  assigned  a  relatively 

large  value.  It  should  be  noted  that  the  two-bond-angle 
term  describes  interactions  between  two  pairs  of  covalently 
bonded  beads  sharing  a  bead  using  a  linear  torsional  spring 
model. 

Intra- sheet  hydrogen-bond  term:  Following  our  prior 
work  on  polyurea  121-23],  the  potential  energy  associated 
with  intra-sheet  inter-chain  hydrogen  bonding  is  expressed 
using  the  18-9  Lennard-Jones  type  of  function  as 
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^hydrogen  _  ^hydrogen 


yiydrogenX  /^ydrogen\  ^ 

-2  ^ 


(3) 


where  is  the  interaction  energy  minimum  and 

represents  the  strength  of  hydrogen-bond  interaction 
between  two  beads,  while  is  the  associated 

hydrogen-bonded  inter-bead  equilibrium  distance. 

Inter-sheet  van  der  Waals  non-bond  interaction  term  As 
far  as  the  non-bond  van  der  Waals  interaction  energy 
involving  beads  residing  on  different  PPTA  sheets  is  con¬ 
cerned,  it  is  represented  using  the  12-6  Lennard- Jones  type 
of  function  as 

„non— bond\  /^non— bond 

(4) 

where  is  the  interaction  energy  minimum  and 

represents  the  strength  of  non-bond  interaction  between 
two  beads,  while  is  the  associated  non-bonded 

inter-bead  equilibrium  distance. 


Coarse-grained  force-field  parameterization 

Procedures  used  for  parameterization  of  each  of  the  four 
parts  of  the  coarse-grained  force-field  are  presented  in  this 
section. 

Intra-chain  bond-stretch  term:  Since  the  coarse-grained 
formulation  used  in  the  present  work  defines  only  one  type 
of  bead,  the  bond-stretch  force-field  term  contains  only  two 
parameters:  and  was  determined  by 

setting  it  equal  to  the  equilibrium  bead/bead  separation,  as 
computed  using  all-atom  molecular  statics  calculations  of 
the  minimum  potential-energy  configuration  of  a  single 
infinitely  long,  perfectly  straight  PPTA  chain. 

To  determine  an  equilibrium  molecular-dynam¬ 

ics-based  free-energy  calculation  procedure  is  combined 
with  the  so-called  ''Adaptive  Biasing  Force''  (ABF)  method 
[22].  During  the  application  of  the  ABF  method,  the  distance 
between  the  centers  of  mass  of  the  (all-atom  represented) 
beads  is  varied  in  small  increments  from  its  initial  (minimum 
energy)  value  by  appropriately  changing  the  relative  position 
of  all  constituent  atoms.  This  was  followed  by  prolonged 
(40-ns  simulation  time,  at  a  temperature  of  300  K)  all-atom 
equilibrium  molecular-dynamics  simulations.  The  trajectory 
results  obtained,  averaged  over  the  simulation  time,  are  used 
to  compute  the  mean  force  acting  on  the  beads,  quantifying 
the  tendency  of  the  beads’  connector  to  recover  its  equilib¬ 
rium  length.  To  cancel  this  force,  a  biasing  force  of  equal 
magnitude  and  opposite  sign  is  equally  partitioned  and 
applied  to  all  the  constituent  atoms  in  the  two  beads.  This 
leaves  the  system  under  zero  mean  force  and  only  the  fluc¬ 
tuating  part  of  the  force  remains,  enabling  the  computation  of 


the  "equilibrium"  thermodynamic  properties  of  the  system, 
e.g.,  mean  potential  energy  at  a  given  value  of  the  bead 
connector  length.  When  the  resulting  potential  energy  versus 
connector  length  results  are  fitted  using  linear  regression 
analysis  in  conjunction  with  Eq  (1),  the  sought-after 
stretching  stiffness  parameter  is  obtained.  The  resulting 
functional  relationship  for  bond-stretch  energy  versus  bon¬ 
ded-bead  separation  is  depicted  in  Fig.  4a. 

Intra-Chain  Two-Bond-Angle  Term  The  two-bond-angle 
force-field  term  also  contains  only  two  parameters: 
and  As  mentioned  earlier,  is  set  to  a  value  of  n. 
As  far  as  the  bending  stiffness  parameter  is  con¬ 

cerned,  it  is  determined  using  an  ABF-based-procedure 
analogous  to  that  described  for  the  case  of  The 

main  difference  is  that  instead  of  a  perfectly  straight  PPTA 
molecule,  a  molecule  consisting  of  two  perfectly  straight 
halves  with  an  included  angle  6  (of  varying  magnitude)  is 
utilized.  The  resulting  functional  relationship  for  bond- 
angle  energy  versus  two-bond-angle  is  depicted  in  Fig.  4b. 

Intra-sheet  hydrogen-bond  term  Due  to  the  presence  of  a 
single  bead  type,  the  hydrogen-bond  part  of  the  coarse¬ 
grained  force-field  is  also  characterized  by  two  parameters: 
^hydrogen  ^ydrogen  two  parameters  Were  deter¬ 

mined  using  the  same  procedure  as  the  ones  described 
above,  except  that  instead  of  a  single  infinitely  long,  per¬ 
fectly  straight  chain,  a  wall  of  parallel,  equally  spaced, 
infinitely  long,  perfectly  straight  chains  is  used,  and  the 
inter-chain  spacing  is  varied.  The  resulting  functional 
relationship  for  hydrogen-bond  energy  versus  intra- sheet 
inter-chain  separation  distance  is  depicted  in  Fig.  4c. 

Inter-sheet  van  der  Waals  non-bond  interaction  term 
The  two  non-bond  interaction  parameters  and 

^non-bond  determined  using  the  same  procedure,  except 
that  an  infinitely  wide/deep/tall  stack  of  PPTA  sheets  (with 
the  ABABAB...  stacking),  each  consisting  of  an  array  of 
parallel,  equally  spaced,  infinitely  long,  perfectly  straight 
chains,  is  used.  The  resulting  functional  relationship  for 
non-bond  energy  versus  inter- sheet  separation  distance  is 
depicted  in  Fig.  4d. 

Parameterization  Summary  A  summary  of  the  eight  coarse¬ 
grained  PPTA  force-field  parameters  is  provided  in  Table  1. 

Coarse-grained  molecular  method 

All-atom  and  coarse-grained  molecular  methods  and  tools 
are  fundamentally  identical,  since  both  treat  materials  as 
assemblies  of  interacting  particles  and  use  force-fields  to 
quantify  various  bonding  and  non-bonding  interactions. 
The  two  main  differences  between  the  two  methods  are 
(a)  the  nature,  and  size/mass,  of  the  constituent  interacting 
(atoms  vs.  beads)  particles  and  (b)  bead/bead  interactions 
are  generally  quite  more  compliant  (softer)  than  their  atom/ 
atom  counterparts. 


^  Springer 


8280 


J  Mater  Sci  (2014)  49:8272-8293 


Fig.  4  Coarse-grained  PPTA  force-field  functions  used  in  the  present  work:  a  bond-stretch  term;  b  two-bond-angle  term;  c  hydrogen-bond  term; 
and  d  non-bond  term 


Table  1  Coarse-grained  PPTA  force-field  parameterization  used  in 
the  present  work 


Parameter 

Unit 

Value 

^^tretch 

kJ/mol  nm^ 

6831.1 

stretch 

nm 

1.31 

^ngle 

kJ/mol  rad^ 

9752 

jangle 

rad 

71 

ghydrogen 

kJ/mol 

22.3 

^hydrogen 

nm 

1.39 

^non-bond 

kJ/mol 

11.1 

non-bond 

^o 

nm 

1.52 

Please  see  Eqs.  (l)-(4)  for  explanation  of  the  symbols 


In  general,  a  complete  formulation  of  a  coarse-grained 
molecular  simulation  problem  includes,  at  a  minimum,  the 
following  four  items:  (a)  a  computational  model  consisting 


of  bead-type  chains  and  larger  structures  (e.g.,  sheets, 
fibrils,  fibers,  etc.);  (b)  a  set  of  inter-bead  force-field 
potentials  which  describe  accurately  various  bonding  and 
non-bonding  bead/bead  interactions;  (c)  a  computational 
method(s)  to  be  used  in  the  simulations;  and  (d)  formula¬ 
tion  of  the  problem  to  be  analyzed.  These  four  items  are 
discussed  in  greater  detail  in  the  remainder  of  this  section. 

Computational  model(s) 

While  the  coarse-grained  computational  approach  allows  a 
substantial  increase  in  the  dimensions  of  the  computational 
cell  analyzed,  a  coarse-grained  analysis  of  a  segment  of  the 
PPTA  fiber  (typical  diameter  on  the  order  of  10  ^m)  would 
still  require  the  use  of  a  prohibitively  large  number  of 
beads.  Considering  the  fact  that  a  majority  of  the  crystal¬ 
line/topological  defects  are  located  in  the  core  of  the  PPTA 
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Fig.  5  Coarse-grained  PPTA  fiber  computational  cell  consisting  of 
nine  coaxial  fibrils  (indicated  using  different  colors):  a  zoomed-out 
view;  and  b  a  close-up  view 

fibers,  and  that  these  defects  play  an  important  role  in  the 
kink-band  nucleation  process  (discussed  later),  the  com¬ 
putational  domain  used  focuses  on  the  core- section  of  the 
PPTA  fibers.  Furthermore,  since  the  transverse  cross-sec¬ 
tion  of  the  computational  domain  is  still  small  in  com¬ 
parison  to  that  of  the  PPTA-fiber  core,  the  computational 
domain  is  assumed  to  extend  indefinitely  in  the  lateral 
directions  (as  well  as  in  the  axial  direction).  The  rectan¬ 
gular-parallelepiped-shaped  computational  cell  used  in  the 
present  work  is  shown  in  Fig.  5a.  The  long  edge,  c,  of  the 
computational  cell  is  aligned  with  the  fiber  axis  (which  is, 
in  turn,  aligned  with  the  z-coordinate  axis).  The  remaining 
two  transverse  edges  (a  and  b)  of  the  cell  are  aligned, 
respectively,  with  the  coordinate  axes  v  and  y.  Typically, 
computational  cells  of  the  size  a  =  24.3  nm,  h  =  24.3  nm, 
and  c  =  28.6  nm  are  used.  The  infinite  extension  of  the 
computational  domain  in  the  three  orthogonal  directions  is 
achieved  by  applying  the  periodic  boundary  conditions 
across  the  computational  cell  faces.  It  should  be  noted  that 
the  dimensions  of  the  computational  cell  used  are  a  com¬ 
promise  between  a  lower  computational  cost  (which  favors 
a  smaller  computational  cell)  and  the  magnitude  of  the 
effect  of  the  computational  cell  size  on  the  material 
strength  (which  favors  a  larger  computational  cell).  In  other 
words,  the  present  computational  cell  sizes  enabled 


analysis  of  the  material  strength/stiffness  behavior  at  a 
relatively  low  computational  cost,  while  ensuring  that 
further  increases  in  the  computational  cell  size  would  not 
yield  significantly  different  quantitative  results. 

Examination  of  the  transverse  cross-section  of  the 
computational  cell  shows  that  the  cell  consists  of  nine 
coaxial  fibrils  (identified  by  different  colors).  While  each 
fibril  contains  the  same  PPTA  crystal  structure,  the  trans¬ 
verse  orientation  of  this  structure  differs  among  the  fibrils. 
However,  in  all  fibrils,  the  orientation  of  the  crystal  lattice 
is  such  that  the  sheet  normal  is  orthogonal  to  the  fiber  axis. 
The  initial  atomic  (or,  more  precisely,  the  bead)  configu¬ 
ration  was  constructed  while  obeying  the  basic  crystallo¬ 
graphic  relationships  within  the  PPTA  phase.  Toward  that 
end,  the  mean  mass  density  within  each  fibril  is  initially  set 
to  be  the  same  (by  properly  adjusting  inter-chain  distances 
within  the  sheets  and  inter-sheet  distances)  for  each  fibril. 
To  equilibrate  the  initial  configuration,  before  conducting 
any  analysis/testing  of  the  computational  cell,  the  potential 
energy  of  the  computational  cell  is  minimized  with  respect 
to  the  spatial  coordinates  of  the  constituent  beads.  It  should 
be  noted  that  the  atomic  configuration  displayed  in  Fig.  5a 
is  the  initial  (before  energy  minimization)  configuration.  A 
close-up  of  this  configuration  is  displayed  in  Fig.  5b.  This 
figure  clearly  shows  the  formation  of  PPTA  sheets  which 
contain  covalent  intra-chain  bonds  and  inter-chain,  intra¬ 
sheet  hydrogen  bonds  (both  depicted  using  the  ball-and- 
stick  representation). 

Force-fields 

The  force-field  functions  used,  the  procedure  employed  for 
their  parameterization  and  the  resulting  parameters’  iden¬ 
tification  were  all  covered  in  “PPTA  coarse-graining  and 
force-field  parameterization”  section  and,  hence,  will  not 
be  repeated  here. 

Computational  method( s ) 

Within  the  coarse-grained  fiber-scale  computational 
investigation  carried  out  in  the  present  work,  two  types  of 
computational  methods,  molecular  statics  (i.e.,  potential- 
energy  minimization)  and  (equilibrium)  molecular 
dynamics  (i.e.,  particle-trajectory  analysis),  are  used.  These 
methods  are  implemented  in  Discover  (a  general-purpose 
atomistic  simulation  program)  from  Accelrys  [24],  the 
computational  tool  which  was  employed  in  the  present 
work.  Within  the  molecular  statics  approach,  the  potential 
energy  of  the  subject  discrete-particle  material  system  is 
minimized  with  respect  to  the  position  of  the  particles/ 
beads.  On  the  other  hand,  within  the  molecular-dynamics 
approach,  gradient  of  the  potential  energy  (as  defined  by 
the  force-field  functions)  with  respect  to  the  particle 
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positions  is  first  used  to  generate  forces  acting  on  the 
particles  and,  then,  the  associated  Newton’s  equations  of 
motion  (for  all  particles)  are  integrated  numerically  in 
order  to  determine  the  temporal  evolution  of  the  particle 
positions  (i.e.,  particle  trajectories).  In  addition,  within  the 
equilibrium  molecular-dynamics  methods,  the  system 
under  consideration  is  coupled  to  an  (external)  environ¬ 
ment  (a  constant  temperature  reservoir,  in  the  present  case) 
which  ensures  that  the  system  remains  in  thermodynamic 
equilibrium  (i.e.,  the  system  does  not  experience  fluxes  of 
thermodynamic  quantities).  In  the  present  work,  an  NVT 
(where  N  is  the  (fixed)  number  of  atoms,  V  is  the  compu¬ 
tational  cell  volume  (also  fixed),  and  T  (=300  K)  is  the 
temperature)  rendition  of  the  molecular-dynamics  method 
is  employed. 

As  will  be  shown  in  the  following,  combined  application 
of  the  molecular  statics  and  equilibrium  molecular- 
dynamics  methods  enables  the  investigation  of  (a)  behavior 
of  the  PPTA  fiber  core  under  axial-compressive/tensile 
loading;  (b)  formation  of  loading-induced  defects  such  as 
kink  bands;  and  (c)  elucidation  of  the  interplay  between 
various  fiber-processing-induced  crystallographic/topolo¬ 
gical  defects  and  the  propensity  of  the  material  to  generate 
additional  defects  (e.g.,  kink  bands  during  loading).  All  the 
calculations  in  the  present  work  were  carried  out  by  cou¬ 
pling  a  user-defined  coarse-grained  force-field  database 
with  the  Discover  solver  [24].  The  procedure  employed  to 
simulate  PPTA-fiber  axial  compression  or  extension  is 
implemented  and  linked  with  Discover  through  the  use  of 
an  input  file  written  in  the  BTCL  script.  This  enabled 
precise  control  of  all  simulation  tasks  such  as  sequential 
contraction  or  expansion  of  the  computational  cell  in  the 
longitudinal  direction  (during  the  simulation  of  the  axial- 
compression/- tension  test). 

Problem  formulation 

As  mentioned  earlier,  the  problem  analyzed  in  the  present 
work  involves  determination  of  the  functional  relationship 
between  the  extent  of  material  degradation  and  long¬ 
itudinal-tensile-strength  loss  in  PPTA  fibers,  as  a  result  of 
the  prior  compressive  or  tensile  axial  loading.  In  a  typical 
coarse-grained  axial  loading  analysis  of  the  PPTA  fibers, 
both  the  preludial  compression/tension  and  the  subsequent 
longitudinal  tension  are  applied  in  small  increments  and 
each  incremental  deformation  is  followed  by  the  compu¬ 
tational  cell  energy  minimization  and  extended  (ca.  40  ns) 
molecular-dynamics  simulation  runs. 

Post-processing  data-reduction  analysis 

The  coarse-grained  fiber-level  particle-trajectory  compu¬ 
tational  results  obtained  are  mainly  used  to  reveal  the 


response  of  single  PPTA  fibers  to  uniaxial  compression/ 
tension  (including  the  microstructural/topological  evolu¬ 
tion  processes  such  as  kink-band  formation  and  fibrilla¬ 
tion).  Fiber  axial  compression/tension  is  simulated  by 
subjecting  the  unit  cell  to  longitudinal  compressive/tensile 
deformation  in  the  z-direction,  while  allowing  for  the 
operation  of  the  Poisson’s  effects  in  the  v-  and  y-directions. 
In  other  words,  at  each  incremental  uniaxial  compression/ 
tension  strain,  the  lateral  dimensions  of  the  cell  are 
adjusted,  while  preserving  the  periodic  boundary  condi¬ 
tions,  until  the  normal  stresses  in  the  lateral  directions 
became  effectively  zero.  In  the  case  of  preludial  axial 
compression/tension,  except  for  using  the  trajectory  data 
for  graphical  display,  no  data  post-processing  procedure 
had  to  be  employed.  However,  the  pre-loaded/unloaded 
computational  cells  were  subsequently  subjected  to  tensile 
reloading  in  order  to  assess  the  extent  of  longitudinal-ten¬ 
sile  strength  (and  Young’s  modulus)  loss  during  preludial 
loading.  In  these  cases  (a)  following  uniaxial  compression/ 
tension  and  the  formation  of  deformation-induced  defects 
such  as  kink  bands  or  fibrillation,  the  unit  cell  is  subjected 
to  longitudinal-tensile  deformation  while  allowing  for  the 
operation  of  the  Poisson’s  effects  in  the  transverse  direc¬ 
tions;  and  (b)  the  trajectory  data  had  to  be  post-processed  to 
extract  and  quantify  the  effect  of  preludial  loading.  During 
tensile  reloading,  the  procedure  described  in  our  recent 
work  [21]  was  utilized  to  compute  the  normal  axial  stress 
and  its  rate  of  change  with  an  increase  in  the  axial  tensile 
strain.  The  longitudinal- tensile  strength  of  the  single  PPTA 
fibril  is  then  set  equal  to  the  value  of  the  corresponding 
normal  stress  at  which  the  rate  of  stress  increases  with 
increase  in  strain  begins  to  decline  ''appreciably''  (denot¬ 
ing  the  onset  of  inelastic-deformation  and/or  damage-ini¬ 
tiation  processes). 

Results  of  the  coarse-grained  analysis 

As  established  in  the  Introduction  section,  PPTA  fibrils/ 
fibers  contain  a  variety  of  synthesis/processing-induced 
crystallographic/morphological  defects  which  control  their 
mechanical  properties.  Consequently,  the  hypothetical  case 
of  crystallographically  perfect  PPTA  fibrils/fibers  has  not 
been  investigated  in  the  present  work.  In  our  recent  work 
[7-11],  a  detailed  overview  was  provided  with  respect  to 
the  character  of  the  main  microstructural/topological 
defects  found  in  PPTA  fibers,  their  classification,  their 
origin  in  relation  to  the  PPTA  synthesis/processing  condi¬ 
tions,  and  their  (estimated)  typical  concentrations.  There¬ 
fore,  a  similar  detailed  overview  will  not  be  provided  here. 
Nevertheless,  it  should  be  noted  that  the  prior  work  [7-11] 
clearly  established  that  among  the  isolated  and  clustered 
point- type  defects,  the  following  has  the  largest  effect  on 
the  PPTA  fibril/fiber  properties  (in  particular,  on  the 
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longitudinal-tensile  strength):  (a)  chain-ends;  (b)  side- 
groups;  and  (c)  voids  and  interstitials.  In  addition,  planar- 
type  defects,  commonly  referred  to  as  defect  bands,  were 
also  found  to  have  a  major  effect  on  the  fibril/fiber  strength. 
In  the  present  work,  only  the  effect  of  (isolated  or  clus¬ 
tered)  chain  ends  or  side-groups  was  investigated.  Chain 
ends  are  the  defects  which  stem  from  the  fact  that  some  of 
the  fibril  chains  do  not  run  continuously  along  the  fiber 
length  but  rather  terminate  (at  one  or  both  ends)  within  the 
fiber.  As  far  as  the  side-group  defects  are  concerned,  they 
are  associated  with  the  presence  of  radicals,  e.g.,  -SO3H, 
-S03“Na^,  between  the  PPTA  sheets.  As  shown  in  our 
previous  work  [7-11],  the  presence  of  these  two  types  of 
defects  and  their  clusters  within  the  (otherwise)  perfect 
PPTA  fibers  causes  local  changes  in  the  material  micro¬ 
structure.  These  micro  structural  changes,  as  will  be  shown 
below,  are  the  root  cause  of  the  reduced  longitudinal-ten¬ 
sile  strength  in  the  as-synthesized  PPTA  fibers. 

Micro  structure  evolution  during  axial  compression 
of  defective  PPTA  fibers 

Evolution  of  the  coarse-grained  micro  structure,  during  axial 
compression,  within  the  computational  cell  containing  (ini¬ 
tially)  as- synthesized  PPTA  fibers  having  chain-end  defects, 
is  depicted  in  Fig.  6a-d.  The  longitudinal  stretch  values 
associated  with  Fig.  6a-d  are  (a)  0.975;  (b)  0.950;  (c)  0.925; 
and  (d)  0.900.  The  lateral  stretches,  in  each  case,  take  on  the 
values  required  by  the  zero-lateral-stress  condition.  Exami¬ 
nation  of  the  results  displayed  in  Fig.  6a-d  clearly  reveals 
formation  and  subsequent  growth/thickening  of  two  counter¬ 
balancing  non-parallel  shear  bands.  In  most  cases,  removal 
of  the  axial-compressive  loading  did  not  result  in  complete 
disappearance  of  the  kink  bands,  indicating  permanent 
damage  induced  into  the  PPTA  fibers  by  the  preludial  lon¬ 
gitudinal  compression.  It  should  be  noted  that  while  the 
results  displayed  in  Fig.  6a-d  appear  to  violate  the  periodic 
boundary  condition  imposed  in  the  lateral  directions,  they 
are  fully  consistent  with  this  type  of  boundary  condition.  In 
other  words,  the  material  protruding  through  the  right  face  of 
the  computational  cell  displayed  in  these  figures  is  missing  at 
the  corresponding  left  section  of  the  computational  cell.  Due 
to  the  chosen  orientation  of  the  unit  cell,  selected  in  such  a 
way  as  to  emphasize  the  formation  of  the  counter-balancing 
shear  bands  protruding  outside  the  right  face  of  the  compu¬ 
tational  cell,  the  corresponding  “void”  portion  of  the  unit 
cell  at  its  left  side  is  not  seen. 

When  isolated  or  clustered  side-group  defects  are 
introduced  into  the  PPTA  fiber,  axial  compression  is  not 
accompanied  by  the  formation  of  well-defined  kink  bands. 
Rather,  micro  structural  changes  resembling  fibrillations 
(flow-induced  chain/fibril  buckling/localized  bending) 
were  observed.  An  example  of  the  microstructural 


evolution  during  axial  compression  of  PPTA  fibers  con¬ 
taining  side-group  defects  is  depicted  in  Fig.  7a-d.  The 
longitudinal  stretch  values  associated  with  Fig.  7a-d  are 
(a)  0.975;  (b)  0.950;  (c)  0.925;  and  (d)  0.900.  The  lateral 
stretches  again,  in  each  case,  take  on  the  values  required  by 
the  zero-lateral-stress  condition.  Examination  of  the  results 
displayed  in  these  figures  reveals  the  formation  of  bulges 
(an  indication  of  the  operation  of  fibrillation  processes)  at 
the  computational-cell  lateral  faces.  The  height  and  the 
breadth  of  these  bulges  increase  in  the  course  of  axial 
compression.  In  most  cases  analyzed,  microstructural 
changes  experienced,  during  preludial  compression,  by  the 
PPTA  fiber  containing  side-group  defects  were  permanent. 
That  is,  with  the  exception  of  minor  elastic -relaxation 
processes,  removal  of  the  compressive  loading  did  not 
restore  the  original/uncompressed  micro  structure. 

Micro  structure  evolution  during  axial  tension  of  defective 
PPTA  fibers 

Evolution  of  the  coarse-grained  micro  structure,  during 
axial  tension,  within  the  computational  cell  containing 
(initially)  as-synthesized  PPTA  fibers  having  chain-end 
defects,  is  depicted  in  Fig.  8a-d.  The  longitudinal  stretch 
values  associated  with  Fig.  8a-d  are  (a)  1.025;  (b)  1.050; 
(c)  1.075;  and  (d)  1.100,  while  the  corresponding  lateral 
stretches  take  on  the  values  required  by  the  zero-lateral- 
stress  condition.  Examination  of  the  results  displayed  in 
Fig.  8a-d  clearly  reveals  that  the  distributions  of  the  lateral 
strains  within  one  of  the  fibrils  are  not  uniform.  Rather, 
there  is  a  localization  of  the  lateral  strains  in  a  region 
adjacent  to  the  pre-existing  chain-end  defect.  In  most  cases, 
removal  of  the  axial-tensile  loading  did  not  result  in 
complete  disappearance  of  the  localized-deformation 
region,  indicating  permanent  damage  induced  in  the  PPTA 
fibers  by  the  preludial  longitudinal  tension. 

When  isolated  or  clustered  side-group  defects  are  pres¬ 
ent  in  the  PPTA  fiber,  axial  tension  does  not  introduce 
observable  non-uniformities  or  defects.  Consequently,  no 
results  pertaining  to  the  microstructural  evolution  during 
axial  tension  of  PPTA  fibers  containing  side-group  defects 
are  shown. 

Longitudinal-tensile  strength  of  the  defective 
as -synthesized  PPTA  fibers 

In  our  prior  work  [7-11],  it  was  clearly  demonstrated  that 
the  longitudinal-tensile  strength  of  the  as- synthesized 
PPTA  fibers  is  highly  sensitive  to  the  presence,  type,  and 
cluster  size  of  the  microstructural/topological  defects.  Due 
to  the  stochastic  nature  of  the  spatial  distribution  of  the 
defect  concentration  within  the  fiber,  the  longitudinal-ten¬ 
sile  strength  of  the  fibers  has  also  been  found  to  be  a 
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(a) 


(c) 


Fig.  6  Evolution  of  the  coarse-grained  microstructure,  during  axial  compression,  within  the  computational  cell  containing  (initially) 
as- synthesized  PPTA  fibers  having  chain-end  defects.  Longitudinal  stretch  values:  a  0.975;  b  0.950;  c  0.925;  and  d  0.900 


stochastic  quantity.  This  finding  is  in  full  agreement  with 
several  experimental  reports  (e.g.,  [25-27]).  In  our  recent 
work  [8],  a  procedure  was  introduced  for  the  derivation  and 
parameterization  of  a  probability  density  function  for  the 
longitudinal-tensile  strength  of  the  as- synthesized  PPTA 
fibers.  Application  of  this  procedure  yielded  the  defective 
as-synthesized  PPTA-fiber  longitudinal- tensile- strength 
probability  density  function,  depicted  in  Fig.  9.  The  mean 
value  of  the  defective  as- synthesized  PPTA  longitudinal- 
tensile  strength,  computed  as  the  first  moment  of  the 
probability  density  function  displayed  in  Fig.  9,  is  found  to 
be  ^5.0  GPa.  This  value  is  in  reasonably  good  agreement 
with  the  measured  PPTA-fiber  longitudinal- tensile  strength 
(typically  found  to  be  in  a  3-5  GPa  range  [e.g.,  23,  26]). 
This  finding  suggests  that  the  procedure  introduced  in  Refs. 
[7-11]  can  reasonably  well  account  for  the  experimental 


observations  regarding  the  PPTA-fiber  longitudinal-tensile 
strength.  Similar  conclusions  can  be  drawn  in  regards  to  the 
elastic  stiffness  moduli  of  the  PPTA  fibers  since  the  pre¬ 
dicted  mean  values  of  the  longitudinal  (345  GPa)  and 
(isotropic)  transverse  (46  GPa)  Young’s  moduli  are  in 
good  agreement  with  their  respective  experimental  coun¬ 
terparts  (350  and  50  GPa)  [8]. 

Residual  longitudinal-tensile  strength  of  defective 
pre-loaded  PPTA  fibers 

The  results  obtained  in  the  present  work  pertaining  to  the 
residual  longitudinal-tensile  strength  of  the  pre-compressed 
or  pre-tensioned  PPTA  fibers  containing  the  two  types  of 
synthesis-induced  micro  structural/morphological  defects 
(i.e.,  chain-end  and  side-group)  established  that  (a)  the 
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(b)  (d) 


Fig.  7  Evolution  of  the  coarse-grained  microstructure,  during  axial  compression,  within  the  computational  cell  containing  (initially) 
as- synthesized  PPTA  fibers  having  side-group  defects.  Longitudinal  stretch  values:  a  0.975;  b  0.950;  c  0.925;  and  d  0.900 


associated  longitudinal-tensile  strength  probability  density 
function  does  not  undergo  a  substantial  change  in  its  shape 
relative  to  its  as- synthesized  counterpart  depicted  in  Fig.  9; 
(b)  the  mean  value  of  the  pre-loaded  longitudinal-tensile 
strength  probability  density  function  is  reduced  relative  to 
its  as-synthesized  counterpart;  and  (c)  the  extent  of  the 
mean  longitudinal-tensile  strength  reduction  is  a  function 
of  both  the  nature  of  the  pre-existing  micro  structural/ 
morphological  defects  (chain-ends  vs.  side-groups)  and  the 
character  of  pre-loading  (compression  vs.  tension). 

Examples  of  the  results  obtained  in  this  portion  of  the 
work  are  shown  in  Fig.  lOa-b.  These  figures  show  the 
percent  loss  of  the  (initial)  mean  fiber  longitudinal-tensile 
strength,  as  a  function  of  the  magnitude  of  the  preludial 
longitudinal  strain.  Figure  10a  pertains  to  the  case  of 
preludial  compression,  while  Fig.  10b  pertains  to  the  case 


of  preludial  tension.  In  both  figures,  two  curves  are  shown: 
one  corresponding  to  the  case  of  defective  fibers  containing 
chain-ends,  while  the  other  corresponds  to  the  case  of 
defective  fibers  containing  side-groups.  It  should  be  noted 
that  non-zero  values  of  the  percent  loss  in  the  fiber  tensile 
strength  are  seen  at  zero  values  of  the  preludial  compres¬ 
sive/tensile  strain.  The  reason  for  this  is  that  the  fabrica¬ 
tion-induced  defects  (i.e.,  chain-ends  and  side-groups) 
present  in  the  PPTA  fiber  give  rise,  themselves,  to  a 
reduction  in  the  material  strength  relative  to  that  found  in 
the  defect-free  material.  In  other  words,  the  values  of  the 
percent  loss  in  the  longitudinal-tensile  strength  are  defined 
with  respect  to  the  longitudinal-tensile  strength  of  defect- 
free  PPTA  fibers.  Examination  and  comparison  of  the 
results  displayed  in  Fig.  lOa-b  reveal  that  (a)  preludial 
longitudinal  compression  gives  rise  to  a  greater  loss  in  the 
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Fig.  8  Evolution  of  the  coarse-grained  microstructure,  during  axial  tension,  within  the  computational  cell  containing  (initially)  as-synthesized 
PPTA  fibers  having  chain-end  defects.  Longitudinal  stretch  values:  a  1.025;  b  1.050;  c  1.075;  and  d  1.100 


fiber  longitudinal-tensile  strength  than  preludial  longitudi¬ 
nal  tension;  (b)  in  the  case  of  preludial  compression,  side- 
groups  appear  to  cause  a  greater  degradation  in  the  fiber 
longitudinal-tensile  strength  than  the  chain-ends;  (c)  con- 
trarily,  in  the  case  of  preludial  tension,  chain-ends  appear 
to  cause  a  greater  degradation  in  the  fiber  longitudinal- 
tensile  strength  than  the  side-groups;  and  (d)  overall,  the 
largest  loss  in  the  fiber  longitudinal-tensile  strength  is 
encountered  in  the  case  of  preludial  longitudinal  com¬ 
pression  and  the  pre-existing  side-group  defects. 

Finite-element  analysis  (FEA)  of  plain-weaving  process 

The  results  obtained  in  the  previous  section  are  used  to 
enrich  the  PPTA-fiber/-yarn  material  model  developed  in 
our  prior  work  [28-30].  Specifically,  the  results  obtained 


are  used  to  parameterize  the  progressive  damage  model 
developed  in  Ref.  [28-30].  Toward  that  end,  the  results 
presented  in  Fig.  lOa-b  are  weighted  by  the  concentrations 
of  the  two  microstructural/topological  defects,  in  order  to 
arrive  at  the  overall  percent  loss  of  longitudinal-tensile 
strength  versus  preludial  (compressive  or  tensile)  strain 
curves.  Fig.  10c.  Following  the  procedure  outlined  in  Refs. 
[28-30],  the  concomitant  loss  in  the  longitudinal  Young’s 
modulus  as  a  function  of  the  preludial  compressive/tensile 
strain  is  next  determined.  The  upgraded  material  model  is 
then  used  in  the  present  finite-element  analysis. 

Problem  formulation 

The  problem  analyzed  in  the  present  work  involves  finite- 
element  modeling  and  analysis  of  a  (simplified)  plain¬ 
weaving  process.  The  objective  of  the  analysis  is  to  assess 
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Fig.  9  Longitudinal-tensile  strength  probability  density  function  of 
the  defective  as- synthesized  PPTA  fibers 


the  extent  of  damage  sustained  by,  and  the  degree  of  lon¬ 
gitudinal-tensile  strength  loss  in,  warp  and  weft  PPTA 
yarns  used  in  this  process  to  manufacture  single-ply 
fabrics. 

Modeling  and  computational  analysis 

The  finite-element  computational  approach  employed  in 
the  present  work  requires  specification  of  the  following: 
(a)  geometrical  model;  (b)  meshed  model;  (c)  computa¬ 
tional  algorithm;  (d)  initial  conditions;  (e)  boundary  con¬ 
ditions;  (f)  steps  of  analysis;  (g)  contact  interactions; 
(h)  material  model(s);  and  (i)  computational  tool.  These 
aspects  are  briefly  overviewed  in  the  remainder  of  this 
subsection. 

Geometrical  model 

Two  views  of  an  example  of  the  geometrical  model/com¬ 
putational  domain  for  the  problem  analyzed  in  this  portion 
of  the  work  are  depicted  in  Fig.  lla-b.  It  is  seen  that  the 
geometrical  model  includes  five  warp  and  five  weft  yams. 
In  the  geometrical  model  shown  in  these  figures,  only  two 
weft  yarns  (labeled  1  and  2)  are  present  in  the  computa¬ 
tional  domain  and  have  been  engaged  by  the  five  warp 
yarns.  In  other  words,  the  warp  yams  are  snaked  about  the 
(pre-tensioned)  weft  yarns  1  and  2,  and  tensioned  at  ±45° 
from  the  fabric  plane  to  create  a  shed  for  the  third  weft 
yarn.  The  remaining  three  weft  yams  are  added  to  the 
model  (one  at  a  time)  after  each  plain- weaving  cycle.  It 
should  be  noted  that  an  oval  cross-sectional  area  of  the 


yarns  was  adopted,  in  the  present  work,  which  is  in 
accordance  with  the  commonly  observed  yarn  cross-sec¬ 
tions  in  the  plain- woven  fabric,  e.g.,  131-34]. 

Meshed  model 

A  close-up  of  the  typical  meshed  model  used  in  the  present 
work  is  depicted  in  Fig.  12.  Each  warp  and  weft  yarn  is 
meshed  using  88560  and  62315  three-dimensional  contin¬ 
uum  eight-node  finite  elements  of  comparable  size, 
respectively.  The  mesh  size  used  in  the  modeling  of  the 
plain-weaving  process  was  determined  by  carrying  out  a 
mesh- sensitivity  analysis  and  represents  a  compromise 
between  the  computational  efficiency  and  accuracy. 

Computational  algorithm 

All  the  calculations  carried  out  in  this  portion  of  the  work 
are  based  on  a  transient,  displacement-based,  purely 
Lagrangian,  conditionally  stable,  explicit  finite-element 
algorithm.  Since  the  plain-weaving  process  is  generally  not 
associated  with  significant  thermal  effects,  such  effects  are 
neglected  in  the  present  work  and,  consequently,  only  an 
isothermal  (ambient  temperature),  purely  mechanical  ana¬ 
lysis  is  carried  out. 

Initial  conditions 

At  the  beginning  of  the  analysis,  all  the  components  of  the 
computational  model  are  assumed  to  be  stationary  and 
stress  free. 

Boundary  conditions 

The  following  boundary  conditions  were  applied  to  the 
computational  domain: 

(a)  Weft  yarns — each  weft  yam  present  in  the  model 
during  a  given  plain- weaving  cycle  is  subjected  to 
the  predetermined  axial  tension.  The  magnitude  of 
this  tension  is  set  to  a  prototypical  value  of  100  MPa 
135];  and 

(b)  Warp  yarns — the  end  face  of  each  warp  yarn 
preceding  the  first  weft  yam  is  kept  fixed,  while 
the  opposite  end  faces  are  subjected  to  the  axial 
tension  of  100  MPa. 

Steps  of  analysis 

The  total  analysis  is  divided  into  a  number  of  discrete 
steps,  each  corresponding  to  a  single  plain-weaving  cycle. 
Within  each  of  these  analysis  steps/cycles,  the  following  is 
carried  out: 
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Fig.  10  Percent  loss  of  mean  fiber  longitudinal-tensile  strength,  as  a 
function  of  the  magnitude  of  the  preludial  longitudinal:  a  compressive; 
and  b  tensile  strain;  and  c  master  curve  relating  the  extent  of  PPTA 

(i)  at  the  beginning  of  each  step,  a  new  weft-yarn  is 
added  to  the  model  and  pre- tensioned; 

(ii)  the  two  families  of  warp  yarns  (inclined  at  +45° 
and  —45°)  are  advanced  in  the  plane  parallel  to 
the  plane  of  the  fabric,  by  one  weft-yarn  spacing; 
and 

(iii)  the  two  warp-yam  families  are  then  swung  in 
order  to  change  the  sign  of  their  inclination 
angle.  This  was  accomplished  by  moving  the  top 
two-pulley  +  belt  mechanism  downward  and  the 
bottom  two-pulley  +  belt  mechanism  upward. 

As  seen  in  Fig.  13,  couplings  are  attached  to  the  ends  of 
(each)  warp  yam.  The  reference  node  of  Coupling  1  has  all 
degrees  of  freedom  fixed,  except  for  rotation  about  the  y-axis 
(the  axis  collinear  with  the  weft  yarns).  The  reference  node 
of  Coupling  2,  on  the  other  hand,  has  two  additional  degrees 


fiber/yam  longitudinal-tensile  strength  loss  as  a  function  of  the 
preludial  compressive/tensile  longitudinal  strain 


of  freedom  (x-  and  z-translations)  free.  Belts  1  and  2  are 
connected  to  the  reference  node  of  Coupling  2.  As  seen  in 
Fig.  1 3 ,  pulling  to  the  right  downward  on  the  free  end  of  Belt 
1  swings  the  respective  family  of  warp  yarns  upward,  while 
pulling  to  the  right  downward  on  the  free  end  of  Belt  2 
swings  the  other  family  of  warp  yarns  downward.  Between 
consecutive  pulling  on  the  two  belts,  a  new  weft  yarn  is 
added  (and  incorporated  into  the  woven  fabric). 

Contact  interactions 

Warp-yam/weft-yam  interactions  are  all  modeled  using  the 
penalty-type  contact  algorithm.  Within  this  algorithm, 
(normal)  penetration  of  the  contacting  surfaces  is  resisted 
by  a  set  of  linear  springs  which  produce  a  contact  pressure 
that  is  proportional  to  the  depth  of  penetration.  Typically, 
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Fig.  11  a  Side  view,  and  b  top  view  of  an  example  of  the  geometrical 
model/computational  domain  for  the  finite-element  analysis  of  the 
plain-weaving  process 


Fig.  12  A  close-up  of  the  typical  meshed  model  used  in  the 
finite-element  analysis  of  the  plain-weaving  process 


maximum  default  values,  which  still  ensure  computational 
stability,  are  assigned  to  the  (penalty)  spring  constants. 
Force  equilibrium  in  a  direction  collinear  with  the  contact 
interface  normal  then  causes  the  penetration  to  acquire  an 
equilibrium  (contact  pressure  dependent)  value.  It  should 
be  noted  that  no  contact  pressures  are  developed  unless 
(and  until)  the  nodes  on  the  ''slave  surface''  contact/pen¬ 
etrate  the  "master  surface."  On  the  other  hand,  the  mag¬ 
nitude  of  the  contact  pressure  that  can  be  developed  is 
unlimited.  As  far  as  the  tangential  contact  interactions 
(responsible  for  transmission  of  the  shear  stresses  across 
the  contact  interface)  are  concerned,  they  are  modeled 
using  a  modified  Coulomb  friction  law.  Within  this  law,  the 
maximum  value  of  the  shear  stresses  that  can  be  trans¬ 
mitted  (before  the  contacting  surfaces  begin  to  slide)  is 


Fig.  13  The  procedure  used  to  carry  out  plain-weaving  of  the  fabric. 
Please  see  text  for  details 


defined  by  a  product  of  the  contact  pressure  and  a  static 
(before  sliding)  and  a  kinetic  (during  sliding)  friction 
coefficient.  In  addition,  to  account  for  the  potential 
occurrence  of  a  "sticking  condition"  (sliding  occurs  by 
shear  fracture  of  the  softer  of  the  two  materials,  rather  than 
by  a  relative  motion  at  the  contact  interface),  a  maximum 
value  of  shear  stress  (equal  to  the  shear  strength  of  the 
softer  material)  that  can  be  transmitted  at  any  level  of  the 
contact  pressure  is  also  specified  [36].  It  should  be  noted 
that,  in  order  to  account  for  the  presence  of  yarn  coating 
(not  modeled  explicitly),  a  low  value  of  the  friction  coef¬ 
ficient  fi  =  0.05  was  used  [8].  It  should  be  also  noted  that 
within  the  yams,  friction  may  develop  between  the  (lightly 
twisted)  contacting  fibers.  Since  the  presence  of  discrete 
fibers  within  individual  yams  is  not  modeled  explicitly,  this 
effect  was  accounted  for  during  the  formulation  of  the 
yarn-level  material  constitutive  model.  This  was  accom¬ 
plished  through  the  use  of  the  so-called  Digital  Element 
Method  and  a  pulled-fiber  test  [8]. 


Material  models 

As  mentioned  earlier,  the  same  material  model  is  used  for 
the  warp-  and  weft-yams.  The  model  was  previously 
developed  in  Refs.  [28-30]  and  upgraded  using  the  coarse¬ 
grained  molecular-level  computational  results  presented  in 
the  previous  section.  The  material  constitutive  model  is  of 
a  linear  orthotropic  elastic  character  with  progressively 
degradable  longitudinal-tensile  strength.  The  coarse-grained 
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molecular-level  computational  results  are  used  to  enrich 
this  material  model  by  making  the  degradation  of  the 
material  tensile  strength  and  Young’s  modulus  be  a  func¬ 
tion  of  the  magnitude  and  the  character  (tension,  com¬ 
pression,  shear)  of  the  prior  and  ongoing  loading.  This 
ensures  that  highly  bent  warp  yarns  at  their  compressive 
side  undergo  substantial  strength/stiffness  degradation 
which  will  prevent  the  development  of  unrealistically  high 
local  stresses.  Since  the  failure  tensile  strain  is  assumed  to 
be  damage  invariant,  the  material  longitudinal  Young’s 
modulus  is  also  treated  as  a  progressively  degradable 
material  property. 

Computational  tool 

The  problem  of  plain-weaving  process  modeling  is  ana¬ 
lyzed  using  an  explicit  solution  algorithm  implemented  in 
ABAQUS/Explicit,  a  general-purpose  finite-element  solver 
[37].  This  algorithm  was  chosen  because  it  is  associated 
with  comparatively  low  computational  cost  when  dealing 
with  three-dimensional  problems  dominated  by  the  contact 
interactions  and  model-change  (i.e.,  progressive  addition  of 
the  weft  yarns).  Since  the  dynamic,  explicit  finite-element 
formulation  is  only  conditionally  stable,  care  had  to  be 
taken  to  ensure  that  the  time  increments  during  the  analysis 
do  not  exceed  the  critical  time  increment  [38]. 

Results 

The  results  presented  in  this  subsection  are  used  to  help 
establish  the  effect  of  differences  in  the  kinematics  and 
deformation  sustained  by  the  warp  and  weft  yarns  during 
plain-weaving  on  their  residual  tensile  strength. 

Material  distribution  during  plain-weaving 

Material  distribution  in  the  course  of  four  consecutive 
cycles  of  the  plain- weaving  process  is  shown  in  Fig.  14a-d. 
For  clarity,  weft  yarns  are  colored  white,  while  warp  yams 
are  colored  red  in  these  figures.  In  Fig.  14a,  it  is  seen  that 
the  first  weft  yam  (labeled  1)  is  completely  engaged,  while 
the  second  (labeled  2)  has  begun  to  be  engaged  by  the 
vertically  moving  two  families  of  warp  yarns.  At  the 
identical  instant  during  the  subsequent  plain-weaving 
cycle.  Fig.  14b,  it  is  seen  that  weft  yarn  2  is  completely 
engaged,  while  weft  yarn  3  has  begun  to  be  engaged  by  the 
vertically  moving  two  families  of  warp  yams.  Figures 
14c-d  show,  respectively,  material  distribution  at  the 
identical  instants  during  the  subsequent  two  weaving  cycles. 
Examination  of  Fig.  14a-d  reveals  that  each  weaving  cycle 
adds  to  the  growing  fabric  section  of  a  length  equal  to  the 
weft-yarn  spacing.  It  should  be  noted  that,  for  improved 
clarity,  the  results  presented  in  Fig.  14a-d  are  generated 


under  relatively  low  tensile  loads  applied  to  the  warp  and 
weft  yarns.  Consequently,  relatively  little  distortion  in  the 
yarn  cross-sectional  area,  except  in  the  warp-yarn/weft-yarn 
contact  region,  is  seen  in  these  figures.  The  extent  of  yarn 
cross-section  distortions  increases  with  the  magnitude  of  the 
tension  forces  applied  to  the  warp  and  weft  yams. 

Weaving-induced  warp -/weft -yarn  longitudinal-tensile 
strength  loss 

Since  the  warp  and  the  weft  yams  are  subjected  to  bending 
and  contact  stresses/strains  during  plain-weaving,  they 
experience  structural  damage  and,  in  turn,  a  loss  in  the 
longitudinal-tensile  strength.  It  should  be  noted  that,  in  the 
present  work,  only  the  normal  contact  stresses/strains  are 
assumed  to  give  rise  to  the  longitudinal-tensile  strength 
loss.  In  other  words,  the  effect  of  tangential  contact  stres¬ 
ses/strains  is  ignored.  This  was  done  for  two  reasons:  (a)  to 
account  for  the  fact  that  fibers/yams  are  generally  coated  to 
minimize  the  friction  coefficient  and,  thus,  sliding  contact 
stresses/strains  and  (b)  to  recognize  the  fact  that  the  coarse¬ 
grained  analysis  carried  out  in  the  present  work  strictly 
modeled  only  the  interior  portion  of  the  PPTA  fibers.  The 
upgraded  material  model  for  the  PPTA  yams,  described 
earlier,  treats  the  material  longitudinal-tensile  strength  as  a 
state  variable,  the  evolution  of  which  is  governed  by  the 
loading  history  of  the  yam.  The  loading  history  of  the  warp 
and  weft  yarns,  in  the  present  case,  is  associated  with  yarn 
crimping  and  contact  accompanying  plain- weaving.  Dis¬ 
tributions  of  the  local  loss  (in  GPa)  in  the  longitudinal- 
tensile  strength  over  a  planar  section  (passing  through  the 
yarn  neutral  axis  and  orthogonal  to  the  fabric  mid-plane)  of 
one  of  the  warp  and  weft  yarns  are  shown  in  Fig.  15a-b, 
respectively.  It  should  be  noted  that  the  yams  depicted  in 
these  figures  reside  in  the  fabric  produced  by  the  plain¬ 
weaving  process  simulated  in  the  present  work.  Examina¬ 
tion  of  the  results  displayed  in  Fig.  15a-b  clearly  reveals 
that  the  warp  yarns  experience  a  larger  extent  of  the  lon¬ 
gitudinal-tensile  strength  loss  and  the  size  of  the  compro¬ 
mised  region  is  larger.  Specifically,  in  the  case  of  the  warp 
yarns,  the  longitudinal-tensile  strength  loss  could  be  as 
high  as  0.28  GPa,  while  in  the  case  of  the  weft  yarns  this 
loss  is  merely  0.08  GPa.  Furthermore,  in  both  types  of 
yarns,  the  maximum  loss  in  the  longitudinal-tensile 
strength  is  observed  in  the  yam  regions  subjected  to 
compression  during  plain- weaving. 

Weaving-induced  warp-/weft-yarn  longitudinal-breaking 
strength  reduction 

In  the  previous  section,  it  was  demonstrated  that  the  warp 
yarns  and,  to  a  lesser  extent  the  weft  yams,  experience  a 
loss  in  their  local  longitudinal- tensile  strength  (modeled  as 
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Fig.  14  Material  distribution  in  the  course  of  four  consecutive  cycles  of  the  plain-weaving  process 


Fig.  16  The  extent  of  failed-element  deletion  prior  to  yarn  failure  in 
the  case  of  the:  a  warp;  and  b  weft  yarns  extracted  from  the  plain- 
woven  fabric  and  subjected  to  uniaxial  tension 


Fig.  15  Distributions  of  the  local  loss  in  the  longitudinal-tensile 
strength  over  a  planar  section  (passing  through  the  yam  neutral  axis 
and  orthogonal  to  the  fabric  mid-plane)  of  one  of  the  a  warp;  and 
b  weft  yarns  residing  in  the  plain-woven  fabric 

a  material  internal  state  variable)  during  plain- weaving.  In 
this  section,  the  yams  depicted  in  Fig.  15a-b  are  extracted 
from  the  fabric  and  subjected  to  axial  loading  (within  a 
computational/virtual  uniaxial-tensile  test)  until  failure. 
This  was  accomplished  by  employing  the  same  type  of 


finite-element  analysis  as  the  one  used  in  the  simulation  of 
the  plain- weaving  process.  Specifically,  couplings  are 
applied  to  the  ends  of  the  extracted  yarns,  and  tensile  forces 
are  applied  along  the  expected  (straight)  axis  of  the  yams. 
To  assess  the  loss  in  the  longitudinal-breaking  strength 
induced  by  weaving,  the  same  virtual  uniaxial-tensile 
testing  is  applied  to  the  as- synthesized  PPTA  yams  [39-42]. 
The  results  obtained  show  that  the  weaving-induced  loss  in 
the  longitudinal-breaking  strength  is  as  follows:  (a)  warp 
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yarns — 10-15  %  and  (b)  weft  yarns — 0-5  %.  These  results 
are  in  reasonable  agreement  with  their  experimental 
counterparts  reported  in  Ref.  [14]  as  (a)  warp  yarns — 
10-25  %  and  (b)  weft  yarns — 0-5  %.  It  should  be  noted 
that  larger  losses  in  the  longitudinal-breaking  strength  than 
in  the  local  longitudinal-tensile  strength  are  observed  in  the 
present  work.  This  is  particularly  true  in  the  case  of  the 
warp  yams.  The  reason  for  this  is  that  the  finite  elements 
within  the  yarn  with  the  most  compromised  longitudinal- 
tensile  strength  become  deleted  in  the  early  stages  of  the 
virtual  uniaxial  testing.  Deletion  of  these  elements  creates 
local  stress-concentration  effects,  leading  to  premature 
failure.  This  is  demonstrated  in  Fig.  16a-b,  in  which  the 
extent  of  failed-element  deletion  prior  to  yam  failure  is 
clearly  greater  in  the  warp-yam  case,  Fig.  16a,  than  in  the 
weft-yarn  case.  Fig.  16b.  It  should  be  further  noted  that,  for 
improved  clarity,  only  short  segments  (containing  the 
region  of  ultimate  failure)  of  the  two  yarns  are  displayed  in 
Fig.  16a-b. 

Conclusions 

Based  on  the  results  obtained  in  the  present  work,  the 
following  summary  remarks  and  main  conclusions  can  be 
drawn: 

1.  The  problem  of  plain- weaving-induced  damage  and 
the  associated  loss  in  the  longitudinal-tensile  strength 
of  the  warp  and  weft  yarns  of  PPTA  is  investigated 
computationally. 

2.  Two  computational  modeling/simulation  techniques  are 
employed:  (a)  coarse-grained  molecular-level  analysis 
is  first  carried  out  to  assess  the  effect  of  preludial  axial  - 
compressive/tensile  (bending)  loading  on  the  longitu¬ 
dinal-tensile  strength  of  the  PPTA  yams  and  (b)  then  a 
conventional  displacement-based  finite-element  ana¬ 
lysis  of  the  plain- weaving  process  is  conducted  in  order 
to  assess  the  extent  of  longitudinal-tensile  strength  of 
the  PPTA  warp  and  weft  yarns  under  common  plain¬ 
weaving  fabric-manufacturing  conditions. 

3.  The  results  obtained  clearly  revealed  that  warp  yams, 
due  to  their  more  extensive  crimping,  sustain  more 
damage  during  the  plain-weaving  process  and,  hence, 
possess  an  inferior  longitudinal-tensile  strength. 

4.  The  computed  results  are  compared  with  the  available 
public-domain  experimental  data  and  a  reasonable 
experiment-to-computation  agreement  is  found. 
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